Bloch oscillations in a periodic lattice under a strong (artificial) 
magnetic field 
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PACS 67 . 85 . -d - Ultracold gases, trapped gases 

PACS 03.65.Vf - Phases: geometric, dynamic or topological 

PACS 73 . 43 . -f - Quantum Hall effects 

Abstract -Inspired by recent developments of artificial gauge fields for atoms and for photons, 
we study the semiclassical dynamics of a quantum particle in a two-dimensional square lattice, 
under the effect of crossed electric and magnetic fields. We provide an elementary derivation of 
how the Berry curvature of the Bloch bands modifies the usual semiclassical equations of motion. 
We then calculate the real- space trajectories of a wave packet and observe that the contribution of 
this anomalous Berry term can be as important as the usual group velocity one. The connection 
with the classical Hall effect is finally clarified. 



In spite of its simple formulation, the coherent dynamics 
of a single quantum particle trapped in a periodic lattice 
potential and subject to a strong magnetic field displays 
a remarkably rich physics and is responsible for several 
important effects of condensed-matter physics. While the 
weak magnetic field case is successfully described by the 
semiclassical theory of Bloch electron as presented in solid- 
state physics textbooks [l] , magnetic field strengths on the 
order of one magnetic flux per unit cell of the lattice are re- 
sponsible for qualitatively novel behaviours. On one hand, 
pioneering work by Hofstadter 2 has anticipated that the 
Bloch bands of the lattice are split into many subbands 
with a fractal structure in the energy vs. magnetic field 
plane. On the other hand, new terms proportional to the 
Berry curvature have to be included in the semiclassical 
theory of transport to account for the peculiar topological 
properties of the Hofstadter bands [3j|4] . 

While many consequences of this physics are nowadays 
known and well understood at a macroscopic level, much 
less experimental literature is available on microscopic 
studies of the basic mechanisms: in ordinary solids, huge 
magnetic fields are required to observe this physics and 
many spurious decoherence phenomena may compete with 
the quantum transport. Only recently, the experimental 
advances in the generation of strong artificial gauge fields 



for (neutral) atoms p\ and for photons p] are opening 
new directions in the study of the quantum dynamics un- 
der strong magnetic fields: first examples in this direction 
are the observation of quantized vortices under the effect 
of a artificial gauge field in an atomic condensate [7] and 
the study of topologically protected photonic edge states 
in photonic crystals (8J. 

In this work we report a theoretical study of Bloch os- 
cillations [9j[l0] of a quantum particle trapped in a two- 
dimensional periodic potential and subject to a strong 
magnetic field orthogonal to the plane of the lattice. Un- 
der the effect of a weak in-plane electric field, the particle 
performs a slow motion across the magnetic Brillouin zone 
that reflects in a peculiar real-space trajectory determined 
not only by the energy dispersion of the band, but also by 
the Berry curvature of the corresponding eigenstates. Re- 
cently, a similar physics was investigated in 11 12] 
atomic lattices of different geometries. In this work 
provide an elementary derivation of the Berry curvature 
terms and we give a complete discussion of the square lat- 
tice case; some interesting features of the connection with 
the classical Hall physics in free space are also clarified. 

The model. — We consider a particle confined in a 
two-dimensional square lattice of spacing a, immersed in 
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a uniform artificial magnetic field B = Bz perpendicular 
to the xy plane of the lattice Q Within the Landau gauge, 
the vector potential associated to the uniform magnetic 
field B has the form A(r) = Bxy. 

Within a standard tight-binding (TB) approximation 
for the periodic lattice potential, the single particle Hamil- 
tonian can be written in the Peierls form 



H 



(id) 



(1) 



as the sum of hopping terms between neighboring sites 
with a non- trivial phase of the hopping amplitude. Here, 
the indices i and j indicate the lattice site, and the sum 
is restricted only to the nearest-neighboring sites. The 
a\ (cLi) operators create (annihilate) one particle at the 
lattice site i and satisfy usual bosonic commutation rules 
[a*, at] = Sij. 

J is a real and positive coefficient quantifying the 
strength of tunneling and the phase of the hopping am- 
plitude from site j to site i can be written in terms of a 
line integral of the vector potential along the hopping path 



A(r) • dr . 



(2) 



While the hopping phase along a single link of the lat- 
tice is manifestly not gauge invariant, the sum of the 
phases accumulated while hopping around a closed con- 
tour is gauge invariant and proportional to the magnetic 
field flux encircled by the contour in units of the flux quan- 
tum </>o = e/27rh. The key parameter determining the 
properties of the model is then a, defined in terms of the 
number of magnetic flux quanta across the elementary cell 
of the lattice as 



^ 1 



2nh 



(3) 



In the Landau gauge used in this work, the hopping am- 



plitude along the x direction has a phase 



0, while 



the one along the y direction has a phase (jx(- = —eBax/h. 

The band structure. — Even though the application 
of a uniform magnetic field does not affect the (discrete) 
translational symmetry of the physical system, this is not 
apparent in the Hamiltonian where the hopping ampli- 
tudes have a spatially dependent phase stemming from the 
spatial dependence of the vector potential A. Of course, 
this breaking of the translational symmetry is only ficti- 
tious, as it is exactly compensated by a suitable gauge 
transformation. 

To diagonalize the Hamiltonian 0, it is useful to con- 
sider the so-called magnetic translation 
whose operators T a are defined as 



group 
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rfi _ — ieAA a (r)-r/h rp 
^a — c ^a 



(4) 



1 As we are dealing with artificial gauge fields, only the products 
eE and eB have a physical meaning. 



in terms of the usual lattice translation operators T a and of 
the vector potential variation AA a (r) = A(r + a) — A(r). 
A complete mathematical characterization of the magnetic 
translation group is given in [l5j[l6]. Here, it is enough 
to note that all T a commute with the Hamiltonian H in 
but they do not in general commute among them. For 
instance, the magnetic translation operators for the two 
unit lattice vectors a. x ^ y in the x, y directions satisfy 



T T — T T 



-2nia 



(5) 



commutation of T ax and T ay is satisfied only if the mag- 
netic flux through a plaquette is an integer multiple of 
the flux quantum <pQ. On the other hand, for a rational 
a = p/q (with p : q co-prime integers), it is possible to de- 
fine a larger magnetic unit cell as the union of q adjacent 
lattice plaquettes: the magnetic flux through the mag- 
netic unit cell is then an integer multiple of 0o, but the 
wider spatial periodicity corresponds to a q times smaller 
magnetic Brillouin zone (MBZ). While the choice of the 
magnetic unit cell is not unique in general, we can take for 
simplicity a magnetic unit cell formed by q adjacent pla- 
quettes along the x direction. With this choice, the MBZ 
becomes of rectangular shape [— ir/qa, 7r/qa]x [— 7r/a, 7r/a], 
with a shorter size along k x . 

Since the Hamiltonian H in ([I]) commutes with the mag- 
netic translation operators, we can look for the common 
eigenstates, which we will call magnetic Block states. To 
conserve the total number of states, the reduced size of the 
magnetic Brillouin zone will be compensated by the sin- 
gle tight-binding band of the non-magnetic a = system 
being folded into q magnetic bands. This was first done 
by Hofstaedter [2Jj who reduced the problem to the solu- 
tion of a one-dimensional difference equation, the so-called 
Harper equation [TF] : 

-J\g(m + l)+g(m-l)] 

— 2J cos(27ram — k y )g(m) = £g(m) . (6) 

Here, the spatial coordinates x and y have been dis- 
cretized on the lattice points as x = ma, y = na, with 
m, n G Z. Thanks to the discrete translational sym- 
metry along the y direction of the vector potential in 
the chosen Landau gauge, a plane wave form can be as- 
sumed for the eigenfunction ip(m,n) along the y direction, 
^(m,n) = ex.-p(ik y n)g(m). 

Even though the eigenvalues of the Hamiltonian ([TJ ar- 
range in a very complex fractal structure in the energy E 
vs. magnetic field parameter a, the so-called Hofstadter 
butterfly, exciting physics is obtained in the case of a ratio- 
nal ct—pjq where the effective potential appearing in the 
Harper equation ([6| is periodic in the x direction with a 
period of q lattice points, in agreement with the magnetic 
translation arguments. 

In particular, for each value of k y , the eigenstates of 
the Harper equation are then classified according to the 
Bloch theorem in terms of a k T wavevector contained in 
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Fig. 1: Energy dispersion of the three magnetic bands in the 
a = 1/3 case in units of the tunneling energy J. 



[—ir/qa,7r/qa]. As a result, we can classify the eigen- 
states of the Hamiltonian ([!]) in terms of k x ^ y within the 
magnetic Brillouin zone and plot the resulting magnetic 
bands. In the non-magnetic a = case, one recovers the 
standard cosine shape of the tight-binding energy band, 
£ = — 2J [cos (A^) + cos {k y )\. In the more interesting 
a = 1/3 case, the three bands have a non-trivial struc- 
ture as shown in Fig. [T] In particular, note the additional 
g-fold periodicity of the Hofstadter bands in the y direc- 
tion. 

For decreasing values of a = 1/g, the standard physics 
of the Landau levels is recovered in that the q bands be- 
come flatter and flatter and the spacing of the lower ones 
tends to the cyclotron frequency uj c = eB/m. This be- 
havior can be recovered from the Harper equation ^ as 
follows: for q — »> oo, the minima of the oscillating potential 
becomes widely separated in space, so tunneling between 
them is suppressed. As the oscillating potential becomes 
smoother, the spacing between the eigenstates decreases 
and the lowest states end up exploring the harmonic bot- 
tom only. More details on this physics will be discussed 
in the last section of the Letter in connection with the 
classical Hall effect. 



Semiclassical dynamics. — In the previous section, 
we have studied the eigenstates of the particle motion in 
the presence of a uniform magnetic field: in particular, we 
have seen that they are classified in terms of a k vector 
within the MBZ. Within a semiclassical picture, we now 
proceed to study the spatial motion of a particle in terms 
of the time evolution of a wave packet contained within 
the n-th magnetic band and localized (as much as the 
Heisenberg principle allows it, of course) in both real and 
k space around r c and k c , respectively. 

According to the modern semiclassical theory of elec- 
tron transport in solids [3j[4] , the time evolution of the 
wavepacket center of mass r c and k c is ruled by the fol- 
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Fig. 2: Normalized Berry curvature Q n /a 2 of the three mag- 
netic bands in the a — 1/3 case. The three panels refer to the 
corresponding panels in Fig[l] We have numerically checked 
that upon integration over the whole MBZ, the Chern num- 
bers for the three bands are respectively C n = —1, 2, —1. 



lowing pair of equations, 

ftk(t) = eE , 
Ht n (t) = V k £ n , k 



eE x n n (k) . 



(7a) 
(7b) 



The first Eq.([7^) describes the usual evolution of the 
(quasi-) momentum under the effect of the electric force 
eE, as discussed in solid-state textbooks [lj. Of course, ac- 
curacy of the semiclassical theory requires that this force is 
weak enough not to induce Landau- Zener interband tran- 
sitions. 

The second Eq.([7j3) is much more interesting: while 
the first term recovers the usual group velocity v^ r k = 
^Vk£ n ,k of a particle living on the energy band £ n ^ the 
second term involves the so-called Berry curvature of the 
band, geometrically introduced as the curvature 



fl n (k) = V k x „4 n ,k 



of the Berry connection 



Ai,k = ^n,k|V k ^n,k) 



(8) 



(9) 



defined in terms of the periodic part u n (k) of the Bloch 
wavefunction VVi,k( r ) = exp(ik • r)i^ n5 k(r). This second 
contribution is the so-called anomalous or Berry velocity, 
which provides a spatial drift proportional to the tempo- 
ral variation of the k vector, and directed orthogonally to 
the applied force. In the solid-state context, this extra 
contribution to the velocity of the electrons was first no- 
ticed by Karplus and Luttinger [|18j in the framework of a 
linear response analysis. More recently, it has been rein- 
terpreted in a more general framework as a Berry phase 
effect in [3|4] . An elementary derivation is sketched in the 
next section. 

Physically, the most striking consequence ofEqs. ([7} is 
that the semiclassical dynamics is not completely encoded 
in the energy spectrum of the Hamiltonian, but keeps track 
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of the relative phase of the eigenstates via their Berry cur- 
vature (J8|. While the Berry curvature identically vanishes 
for systems with simultaneous time-reversal and spatial 
inversion symmetries, a non-trivial value of H n ,k stems in 
our case from the presence of the external magnetic field 
breaking time-reversal. 

A plot of the Berry curvature for the three magnetic 
bands in the a = 1/3 case is shown in Figj2j its general 
structure reflects the g-fold degeneracy of the energy bands 
in the k y direction. Thanks to the inversion symmetry of 
the system, the curvature is symmetric under k — » — k. 
Finally, the sum of the Berry curvature of the different 
bands at the same point k in momentum space identically 
vanishes throughout the whole MBZ. 

The most celebrated consequence of the Berry curva- 
ture in electron transport is the Hall effect, in particular 
the so-called integer quantum Hall effect: the quantized 
value of the Hall conductivity can be related to the Chern 
number of the occupied bands, an integer- valued topologi- 
cal invariant defined as the integral of the Berry curvature 
over the magnetic Brillouin zone in units of 2tt 19 . Other 



electronic systems displaying bands with non-trivial Berry 
curvature are ferromagnetic Fe crystals 20 or graphene 



single layers with a staggered potential breaking the spa- 
tial inversion 
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Elementary derivation of the anomalous Berry 
velocity. — Before proceeding, it is worth presenting an 
elementary derivation of the semiclassical equations of mo- 
tions that does not rely on the Lagrangian formalism used 
in the original works 3,4. For simplicity, we restrict to 
the simplest case of a weak and constant electric field E. 
With a suitable choice of gauge, both the magnetic and 
the electric fields can be included in a (slowly varying) 
vector potential 



A(r,t) = A B (r)-Et, 



(10) 



the first constant term Ae(r) referring to the constant 
magnetic field, V x Ab = B. The system Hamiltonian 
H(t) will have again the same form ([!]) but the hopping 
phases defined by ([2| in terms of the vector potential ( [To] ) 
will be now time-dependent. The Hamiltonian S(t) has 
the same magnetic translation symmetry properties as the 
Hamiltonian and the eigenstates can be classified in 
terms of the momentum k within the MBZ. Provided the 
temporal variation A(t) is slow enough, the adiabatic the- 
ory guarantees that no inter-band transition will occur. 
Thanks to the magnetic translational symmetry of H, each 
k component will be decoupled from all others and will 
just acquire a k-dependent phase factor. 

Under this approximation, the semiclassical state will 
be written in the form 



\w(t)) = / d 2 k™(k)|?A 



(ii) 



where k) ^ s ^ ne e ig ens tate of the Hamiltonian H(t) at 
a wavevector k on the magnetic band n of interest at time 



t. During time-evolution, the eigenfunctions of the Hamil- 
tonian acquires a k-dependent phase 0(k, t), according to 
Berry's adiabatic theorem, satisfying 



d 



R^(k,t) = -^ + eE-4 t , ) k: 



(12) 



in terms of the eigenenergy S n \ of H (t) at the given time 
t. The physical meaning of the Berry term 



(*) 

n,k 



dt 



(13) 



becomes transparent once we note that the eigenstate of 
the time-dependent S(t) at wavevector k coincides with 
the one at wavevector k + eEt of the original Hamilto- 



^k+eEt and 



nian H = H(0) as defined in |l|, £ q 

IV^ n ,k) = l^n k+eEt)- ^his can ^ e easn y shown by means 
of a gauge transformation to eliminate the time-dependent 
component of the vector potential. Analogously, the Berry 
term is related to the usual Berry connection on the MBZ 

by Ai,k = An, k+eEt- 



ky A 



k+ dk , k+hk+q (t) k+bk+q(t+dt) 



k+q(t) k+q(t+3t) 



Fig. 3: Sketch of the integration paths involved in the time 
evolution of the phase 0(k, £). We indicate q(t) = eEt. 

To obtain the evolution equation for the real space po- 
sition r c of the wavepacket, it is then enough to apply the 



stationary-phase approximation to the wavefunction (11) 



using the form (12) of the phase. Under this approxima- 



tion, the position r c is determined by the k-gradient of the 
phase 0(k, t) as 



- V k [ArgHk)] + 0(k, t) + AnM+em] 



(14) 



note in particular the non-trivial term proportional to the 
Berry connection A that accounts for the k-dependence of 
the phase of the Bloch eigenstate. The equation for the 
velocity dr c /dt is then obtained as a time derivative 



dr c 
dt 



:E X Q nM , 



(15) 



the first term is a direct consequence of the definition 
v nk = H ^k£n,k- The term involving the Berry curvature 
requires a bit more care and can be intuitively understood 
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from the sketch in FigJ3J the difference of the phases ac- 
quired after a time t by the components of wavevectors k 
and k + 5k is given by the difference of the integrals of the 
Berry connection A taken along the two horizontal sides 
of the parallelogram. The variation of the Berry connec- 



tion term in (14) corresponds instead to the difference of 



the integral of the Berry connection taken along the two 
shorter sides. The sum of all these contributions recovers 
the time-derivative of the integral of the Berry connec- 
tion around the whole parallelogram, that is the flux of 
the Berry curvature through the infinitesimal gray-shaded 
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[E5t x 5k] n = -(E x St) 5k. 



(16) 



Magnetic Bloch oscillations. — A most fascinating 
feature of the quantum dynamics of a particle in a periodic 
lattice are the so-called Bloch oscillations (BO), that is a 
periodic motion in both real and momentum space under 
the effect of a uniform and time-independent force. In the 
solid-state context, scattering and dephasing effects due to 
impurities in the crystal and to interaction with phonons 
make the observation of BO's in ordinary solids extremely 
difficult unless artificial superlattices with enhanced nm 
periodicities are used 9 . On the other hand, the ex- 
tremely clean periodic potential experienced by atoms in 
optical lattices has allowed for a neat observations of BO's, 
mostly in momentum space 1 10 22|[23] . Very recently, this 
same effect has been exploited to map the band structure 
of atoms in a honeycomb-like lattice potential [24] . 

Starting from the semiclassical Eq.([7^), the oscillation 
frequency along each principal direction x, y of the lattice 
is proportional to the corresponding component eE x y of 



eE x ^ y a/27rh: the momentum 



the applied force, = 
k wraps around the FBZ at a constant speed set by the 
applied force; whenever it reaches an edge of the FBZ, a 
Bragg scattering process makes it to reappear on the other 
side. 

Under the same semiclassical approximation, the spa- 
tial motion is ruled by Eq.([7j3): in standard lattices with 
zero Berry curvature ft = 0, the shape of the BO is fully 
determined by the group velocity term. In the simplest 
case where the force is directed along a principal axis and 
the momentum starts from a high symmetry point like 
k = 0, the BO have the usual back-and-forth shape along 
the direction of the force and their spatial amplitude is 
Ax = A£/\eE\ where A£ is the width of the energy band 
explored by the particle momentum along its path across 
the FBZ. In the general case, the complicate shape of the 
band energy £ n ^ across the FBZ leads to complex trajec- 
tories in real space, which do not typically close on them- 
selves and can show lateral drifts 
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The situation is far more interesting in the presence of a 
finite magnetic field a > 0. To concentrate on the physics 
of interest, we again start from the high symmetry k = 
point in k-space and we direct the applied force along 
the symmetry direction x\ as a result, the wave vector 



15 
10 
5 

-5 
-10 \ 
-15 



MBO 1st band 

MBO 2nd band 

MBO 3rd band 




-50 -40 -30 -20 



-10 

x/a 



10 20 30 



Fig. 4: Trajectory of the magnetic Bloch oscillations (MBO) in 
the a = 1/3 case. The applied electric field is in the x direction 
and has a strength 0.01 J /a. The starting point of the MBO is 
at the origin of both real and reciprocal spaces, r = 0, k = 0. 



k sweeps the MBZ along the k y = line and the group 
velocity remains by symmetry always directed along x. 
Still, as it is illustrated in Figj4]the particle motion from 
the initial position r = shows a non- vanishing drift along 
the y direction due to the Berry curvature in addition to 
the periodic Bloch oscillation along x: the different curves 
show these magnetic Bloch oscillations (MBO) for each of 
the three magnetic bands of the a = 1/3 case, whose Berry 
curvatures differ in both magnitude and sign. The relative 
amplitude of the motion along x and y is determined by 
the strength of the applied force E: the stronger this force, 
the shorter the amplitude of the periodic motion along x 
and the more important the lateral drift term in the y 
direction. 

In the general case of an arbitrary initial wavevector k, 
the trajectory of the MBO is much more complicate, which 
makes it difficult to isolate the contribution of the Berry 
curvature. A clever time-reversal protocol to overcome 
this difficulty and obtain a complete mapping of the Berry 
curvature over the whole MBZ was proposed in [12] . 

Low magnetic field regime and Hall effect. — In 

order to make our discussion complete, it is interesting 
to see how the classical Hall effect for a charged particle 
in free space under crossed electric and magnetic fields 
is recovered in our formalism. In order to recover this 
physics in the lattice model, it is enough to let the lat- 
tice spacing a — » at constant values of the effective mass 
m* = h 2 /2Ja 2 and of the applied electric E and mag- 
netic B fields. In this limit, the discreteness of the lat- 
tice should disappear and the system should recover the 
textbook classical mechanics prediction: the particle mo- 
tion follows a cycloidal trajectory resulting from the su- 
perposition of a circular motion at the cyclotron frequency 
uj c = eB/m and a uniform drift at a Hall velocity 



w H = -E x B/B 2 . 



(17) 
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In a most important limiting case, the circular motion de- 
generates into a point and one is left with only the uniform 
drift at the Hall velocity v#. 



3 




the different magnetic bands, 
the cyclotron frequency uj c . 



that are indeed spaced by 



Conclusions. — We have studied the semiclassical 
dynamics of a quantum mechanical particle in a two- 
dimensional lattice in the presence of a strong (artificial) 
magnetic field. Under a constant (artificial) electric field, 
the main effect of the Berry curvature of the magnetic 
bands is to introduce a lateral drift in addition to the 
periodic Bloch oscillations along the electric field direc- 
tion. An explicit derivation of this result is presented in 
terms of elementary quantum mechanics. In the contin- 
uum limit, the amplitude of the Bloch oscillations tends 
to zero, while the lateral drift due to the Berry curvature 
recovers the classical Hall velocity in crossed electric and 
magnetic fields. 



Fig. 5: Real-space trajectories (measured in units of the mag- 
netic length Im = y/h/eB) of a particle under crossed electric 
and magnetic fields for a magnetic field strength a = 1/9. The 
black solid line is the classical trajectory in free space for an 
initial condition such that the circular component of the mo- 
tion degenerates to a single point. The red dashed line is the 
MBO in the lowest magnetic band. Corresponding symbols on 
top of each trajectory indicate the positions at corresponding 
times. The blue dashed line is the BO trajectory for vanishing 
magnetic field a — 0. The inset is a magnified view of the 
region close to x = 0. 

In Figj5]we show how this behavior is recovered in our 
quantum mechanical lattice model. For low values of the 
magnetic field strength a, the lowest magnetic energy 
bands typically flatten tending to a series of equispaced 
Landau levels separated by the cyclotron frequency u c in 
the a — » limit. Correspondingly, also the Berry cur- 
vature °f the lowest band tends to a constant across 
the whole FBZ. Approximating Qq with a constant and 
knowing that the Chern number of this band is Co = — 1, 
one can heuristically determine its value to ffco = —hz/eB. 
Inserting these remarks into the second semiclassical equa- 
tion Eq.([7]3), one immediately sees that the flatness of the 
bands reflects into a very small amplitude of the MBO, 
while the Berry curvature term recovers the drift at the 
Hall velocity v#: it is quite remarkable to see that to 
recover the Hall effect in the quantum mechanical lattice 
model one has to properly keep into account the Berry cur- 
vature! The small-sized MBO that are visible in the inset 
are due to the exponentially small, but still finite width 
of the magnetic bands: it is important to note that their 
frequency is determined by the strength E of the applied 
electric field and are distinct from the cyclotron motion 
of classical mechanics, whose frequency is determined by 
the strength B of the magnetic field. For this latter to be 
recovered in the quantum picture, one should in fact take 
as the initial state of the particle a linear superposition of 



We thank R. O. Umucalilar for useful discussions. 
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